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O 

f**^ , Abstract. The purpose of this article is to give a streamlined and sclf- 

^\1 ■ contained treatment of the long-time asymptotics of the Toda lattice for de- 

' caying initial data in the soliton and in the similarity region via the method 

j^ ' of nonlinear steepest descent. 






1. Introduction 
The simplest model of a solid is a chain of particles with nearest neighbor inter- 



(1.1) n{p, q) = Yl (^^^ + V{q{n +l,t)- q{n, t))) , 



(N 

f^ ■ action. The Hamiltonian of such a system is given by 

'p{n,ty 

nez 

^ , where q{n, t) is the displacement of the n-th particle from its equilibrium position, 

(y^ ' p[n,t) is its momentum (mass m = 1), and V{r) is the interaction potential. 

C^ . Restricting the attention to finitely many particles (e.g., by imposing periodic 

boundary conditions) and to the harmonic interaction V{r) = ^, the equations of 
motion form a linear system of differential equations with constant coefficients. The 

'bl \ solution is then given by a superposition of the associated normal modes. Around 

1950 it was generally believed that a generic nonlinear perturbation would yield 

f~^ \ to thermalization. That is, for any initial condition the energy should eventually 

be equally distributed over all normal modes. In 1955 Enrico Fermi, John Pasta, 
and Stanislaw Ulam carried out a seemingly innocent computer experiment at Los 
Alamos, |15| . to investigate the rate of approach to the equipartition of energy. 



> 

X. 

H , However, much to everybody's surprise, the experiment indicated, instead of the 

- expected thermalization, a quasi-periodic motion of the system! Many attempts 

were made to explain this result but it was not until ten years later that Martin 
Kruskal and Norman Zabusky, [33]) revealed the connections with solitons (see [5] 
for further historical information and a pedagogical discussion). 

This had a big impact on soliton mathematics and led to an explosive growth 
in the last decades. In particular, it led to the search for a potential V{r) for 
which the above system has soliton solutions. By considering addition formulas for 
elliptic functions, Morikazu Toda came up with the choice V{r) = e^*" + r — 1. The 
corresponding system is now known as the Toda equation, j40| . 
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The equation of motion in this case reads exphcitly 



d dH{p,q) 

—p[n,t) = -— — — — 
at oq[n,t) 



^ Q--iqin.t)-qin-l.t)) 



^-(g(n+l,t)-g(«,t)) 



The important property of the Toda equation is the existence of so called soliton 
solutions, that is, pulslike waves which spread in time without changing their size 
or shape and interact with each other in a particle-like way. This is a surprising 
phenomenon, since for a generic linear equation one would expect spreading of waves 
(dispersion) and for a generic nonlinear force one would expect that solutions only 
exist for a finite time (breaking of waves). Obviously our particular force is such 
that both phenomena cancel each other giving rise to a stable wave existing for all 
time! 

In fact, in the simplest case of one soliton, you can easily verify that this solution 
is given by 

1 + iTT^^^TT exp(— 2Kn + 2(7sinh(K)t) \ 



(1.3) qi{n,t) =q++\og 



1 



■ exp(— 2K(n + 1) + 2crsinh(K)t) 



with K,7 > and a G {±1}- It describes a single bump traveling through the 



Figure 1. One soliton qi{n,0) with k = 1, 7 = 1, and qo — 0. 



crystal with speed cr sinh(K)/K and width proportional to 1/k. In other words, the 
smaller the soliton the faster it propagates. It results in a total displacement 2k of 
the crystal. 

However, this is just the tip of the iceberg and can be generalized to the 7V-soliton 
solution 



(1.4) 

where 

(1.5) CN{n,t)^ 



/ det(I + C^v(M)) A 
q,in,t) = ,^ +log ( det(I + C.(n + M)) j ' 



V7*(Mhj(M) 



^j{n,t) =7jC 



-2K,jn — 2(7j sinh{Kj)t 



\ / l<i,j<N 

with Kj,jj > and aj £ {±1}- The case A^ = 1 coincides with the one soliton 
solution from above and asymptotically, as i — ^ 00, the A^-soliton solution can be 
written as a sum of one-soliton solutions. 

Historically such solitary waves were first observed by the naval architect John 
Scott Russel [3S], who followed the bow wave of a barge which moved along a 
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Figure 2. Numerically computed solution q{n, 150) of the Toda 
lattice, with initial condition all particles at rest in their equilib- 
rium positions except for the one in the middle which is displaced 
byl. 

channel maintaining its speed and size (see the review article |33j for further infor- 
mation) . 

The importance of these solitary waves is that they constitute the stable part 
of the solutions arising from arbitrary short range initial conditions and can be 
used to explain the quasi-periodic behaviour found by Fermi, Pasta, and Ulam. In 
fact, the classical result discovered by Zabusky and Kruskal [43] states that every 
"short range" initial condition eventually splits into a number of stable solitons 
and a decaying background radiation component. This is illustrated in Figure [2] 
which shows the numerically computed solution q{n, t) corresponding to the initial 
condition g(n,0) = (5o,ra, p{n,Q) = at some large time t = 130. You can see the 
soliton region | j | > 1 with two single soliton on the very left respectively right and 
the similarity region | j | < 1 where there is a continuous displacement plus some 
small oscillations which decay like t~^''^ and are asymptotically given by 

(1.6) q{n,t)-2\og{n{zo))+(--^^^^) cos(t$o(;2o) + K^o)log(t)-<5(zo)), 

where zq — e'^" is a slow variable depending only on j and the functions Tq{zo), 
iy{zQ), $0(^0), and S{zo) are explicitly given in terms of the scattering data asso- 
ciated with the initial data. Our main goal will be to mathematically justify this 
formula for the solution in the similarity region | j | < 1 (Theorem 12. 2p and to 
show that the solution splits into a number of solitons in the soliton region | j | > 1 
(Theorem [2H). 

Existence of soliton solutions is usually connected to complete integrability of 
the system, and this is also true for the Toda equation. To see that the Toda 
equation is indeed integrable we introduce Flaschka's variables (TB] 

(1.7) a{n,t) = lc-(i{n+i,t)-q{n,t))/2^ ^(^^^) ^ --p{n,t) 
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and obtain the form most convenient for us 

f^ait)^a{t){b+{t)-bit)), 

(1.8) j^b{t)=2[aitf-a-itr). 

Here we have used the abbreviation 

(1.9) /±(n)-/(n±l). 

Note that if q{n, t) ~> q± sufficiently fast as n -^ ±oo, the converse map is given 

by 

(1.10) q{n,t) = q+ + 2log I f[{2a{j,t))\, p{n,t) ^ -2b{n,t). 

Moreover, q{n,t) — > q±, p{n,t) — > as |n| — > cx) corresponds to a{n,t) — > ^ 
b{n,t) ^0. 

To show complete integrability it suffices to find a so-called Lax pair [27 , that 
is, two operators H{t), P{t) in £^(Z) such that the Lax equation 

(1.11) ^^H{t)^P{t)H{t)-H{t)P{t) 



is equivalent to (|1.8p . One can easily convince oneself that the right choice is 
H{t) = a{t)S+ + a-{t)S- + 5(f), 

(1.12) P{t) = a{t)S+ ~a-{t)S-, 

where {S^f){n) = f^{n) = J{n± 1) are the shift operators. Now the Lax equation 
(jl.lip implies that the operators H{t) for different t € M are unitarily equivalent 
(cf. [Ml Thm. 12.4]): 

Theorem 1.1. Let P(t) be a family of bounded skew-adjoint operators, such that 
t H-> P{t) is differentiate. Then there exists a family of unitary propagators U(t, s) 
for P(t), that is, 

(1.13) ^U(t,s)=P(t)U(t,s), U(s,s)=I. 
dt 



Moreover, the Lax equation (ll.ll[) implies 

(1.14) H{t) = U{t, s)H{s)U{t, s)-\ 

This result has several important consequences. First of all it implies global 
existence of solutions of the Toda lattice. In fact, considering the Banach space of 
all bounded real-valued coefficients (a(n), b{n)) (with the sup norm), local existence 
follows from standard results for differential equations in Banach spaces. Moreover, 
Theorem 11.11 implies that the norm ||_ff(i)|| is constant, which in turn provides a 
uniform bound on the coefficients of H{t), 

(1.15) ||a(i)|U + IIKOIloo < 2\\H{t)\\ = 2\\H{0)\\. 

Hence solutions of the Toda lattice cannot blow up and are global in time (see [351 
Sect. 12.2] for details). 



LONG-TIME ASYMPTOTICS OF THE TODA LATTICE 5 

Second, it provides an infinite sequence of conservation laws expected from a 
completely integrable system. Indeed, if the Lax equation (ll.ll|l holds for H{t), it 
automatically also holds for H{ty . Taking traces shows that 

(1.16) u{H{ty~Hl), jGN, 

is an infinite sequence of conserved quantities, where Hq is the operator correspond- 
ing to the constant solution ao(7i, t) = ^j ^^{n, t) = (it is needed to make the trace 
converge). Introducing a suitable symplectic structure, they can be shown to be in 
involution as well ([IHl Sect. 1.7]). For example, 

tr {H{t) -Ho)=J2 ^("' *) = - ^ II P("' *) ai^d 

(1.17) tr {H{tf - Hi) = J2 Hn, tf + 2{a{n, tf - 1) = ]^n{p, q) 

nez 
correspond to conservation of the total momentum and the total energy, respec- 
tively. 

These observations pave the way for a solution of the Toda equation via the 
inverse scattering transform originally invented by Gardner, Green, Kruskal, and 
Miura [17] for the Korteweg-De Vries equation (see '38j Sect. 13.4] for the case 
of the Toda lattice). In particular. Theorem 11.11 implies that the operators -ff(i), 
i g M, are unitarily equivalent and that the spectrum a{H{t)) is independent of 
t. Now the general idea is to find suitable spectral data S{H{t)) for H{t) which 
uniquely determine H{t). Then equation p. lip can be used to derive linear evolu- 
tion equations for S{H{t)) which are easy to solve. In our case these data will be 
the so called scattering data and the formal procedure (which can be thought of as 
a nonlinear Fourier transform) is summarized below: 



time evolution 
S{H{0)) S{H{t)) 



direct 
scattering 



inverse 
scatterina 



(a(0),6(0)) {a{t),b{t)) 

The inverse scattering step will be done by reformulating the problem as a Riemann- 
Hilbcrt factorization problem. This Riemann-Hilbert problem will then be analyzed 
using the method of nonlinear steepest descent by Deift and Zhou [4] (which is the 
nonlinear analog of the steepest descent for Fourier type integrals). In fact, one of 
our goals is to give a complete and expository introduction to this method. We are 
trying to present a streamlined and simplified approach with complete proofs. In 
particular, we have added two appendices which show how to solve the localized 
Riemann-Hilbert problem on a small cross via parabolic cylinder functions and how 
to rewrite Riemann-Hilbert problems as singular integral equations. Only some ba- 
sic knowledge on Riemann-Hilbert problems, which can be found for example in 
the beautiful lecture notes by Deift [3], is required. 
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For further information on the history of the steepest descent method, which was 
inspired by earher work of Manakov |28, and Its 19J, and the problem of finding 
the long-time asymptotics for integrable nonlinear wave equations, we refer to the 
survey by Deift, Its, and Zhou [7j. 

More information on the Toda lattice can be found in the monographs by Faddeev 
and Takhtajan [lA,, Gesztesy, Holden, Michor, and Teschl [18], Teschl [38|, or Toda 
[40j . Here we partly followed the review article [39j . A much more comprehensive 
guide to the literature can be found in Section 1.8 of [TSj . 

First results on the long-time asymptotics of the doubly infinite Toda lattice 
were given by Novokshenov and HabibuUin 32J and Kamvissis [20 . Long-time 
asymptotics for the finite and semi infinite Toda lattice can be found in Moser [30] 
and Deift, Li, and Tomei ^, respectively. The long-time behaviour of Toda shock 
problem was investigated by Kamvissis [21] and Venakides, Deift, and Oba ^41j and 
of the Toda rarefaction problem by Deift, Kamvissis, Kriecherbauer, and Zhou [llj . 
For the case of a periodic driving force see Deift, Kriecherbauer, and Venakides [8]. 

Finally, we also want to mention that one could replace the constant background 
solution by a periodic one. However, this case exhibits a much different behaviour, 
as was pointed out by Kamvissis and Teschl in [22] (see also [12], [13], [23], [24], 
and [26_ for a rigorous mathematical treatment). 

2. Main results 

As stated in the introduction, we want to compute the long-time asymptotics 
for the doubly infinite Toda lattice which reads in Flaschka's variables 

b{n,t) = 2{a{n,t)^-a{n-l,tf), 
a{n, t) — a{n, t){b{n + 1, i) — b{n, i)), 

(n, i) e Z X R. Here the dot denotes differentiation with respect to time. We will 
consider solutions (a, b) satisfying 

(2.2) ^(1 + \n\y+\\a{n,t) - ^\ + \b{n,t)\) < cx) 

71 

for some ? G N for one (and hence for all, see [35]) i G M. It is well-known that 
the corresponding initial value problem has unique global solutions which can be 
computed via the inverse scattering transform '38]- 

The long-time asymptotics were first derived by Novokshenov and HabibuUin [32] 
and were later made rigorous by Kamvissis [5D] under the additional assumption 
that no solitons are present. The case of solitons was recently investigated by us in 

m- 

As one of our main simplifications in contradistinction to [20) we will work with 
the vector Riemann-Hilbert problem which arises naturally from the inverse scatter- 
ing theory, thus avoiding the detour over the associated matrix Riemann-Hilbert 
problem. This also avoids the singularities appearing in the matrix Riemann- 
Hilbert problem in case the reflection coefficient is —1 at the band edges. 

To state the main results, we begin by recalling that the sequences a(n, f), &(n, i), 
n e Z, for fixed f e M, are uniquely determined by its scattering data, that is, by 
its right reflection coefficient R+{z, t), \z\ — 1, and its eigenvalues Xj <E (— oo, —1) U 
(l,oo), j = l,...,iV, together with the corresponding right norming constants 
7+j(t) > 0, j = 1, . . . , A^. It is well-known that under the assumption (|2.2p the 
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reflection eoeflicients are C'+^(T). Rather than in the complex plane, we will work 
on the unit disc using the usual Joukowski transformation 



(2.3) 



2 I z 



z^X- v/A2 - 1, 



Ae 



UKl. 



In these new coordinates the eigenvalues A^ G (— oo, — 1) U (1, oo) will be denoted by 
Q G (—1,0) U (0, 1). The continuous spectrum [—1, 1] is mapped to the unit circle 
T. Moreover, the phase of the associated Riemann-Hilbert problem is given by 



2j\og{z) 



(2.4) $(z) = z-z" 

and the stationary phase points, $'(z) = 0, are denoted by 



(2.5) 



zo 



n 
1 



i^r 



n 
1 



<?'' 



and correspond to 
(2.6) 



n 
1' 



Here the branch of the square root is chosen such that Iva^-s/z) > 0. For j < — 1 
we have zq € (0,1), for — 1 < j < 1 we have za € T (and hence Zq^ — zq), 
and for J > 1 we have zq £ (—1,0). For \j\ > 1 we will also need the value 
Co e (-1, 0) U (0, 1) defined via Rc($(Co)) = 0, that is. 



(2.7) 



n 
1 



Co - c, 

21og(|Co|)' 



We will set Co = ^1 if If I ^ 1 for notational convenience. A simple analysis shows 
that for J < — 1 we have < Co < 2^0 < 1 and for f > 1 we have — 1 < 2:0 < Co < 0- 
Furthermore, recall that the transmission coefficient T{z), \z\ < 1, is time inde- 
pendent and can be reconstructed using the Poisson- Jensen formula. In particular, 
we define the partial transmission coefficient with respect to zq by 



T{Z, Zq) 



(2.8) <^ 



n lai^, 

ae(Co,o) 

n |ai^)exp('^7log(|r(.)|)f±|^ 



^oe(-i,o), 

\zo\ - 1, 



n lai^ expU./log(|T(s)|)f±ft , zoe(0,l). 

I, Vae(-i,o)u(Co,i) ' / \ T / 

Here, in the case zq € T, the integral is to be taken along the arc T,{zo) ~ {z £ 
T| Re(z) < Re(zo)} oriented counterclockwise. For zq £ (—1,0) we set ^(^o) — 
and for zq £ (0, 1) we set S(zo) = T. Then T{z, zq) is meromorphic for z £ C\E(zo). 
Observe that T{z,zo) = T{z) once zo G (0, 1) and (0,Co) contains no eigenvalues. 
Moreover, r(z, zo) can be computed in terms of the scattering data since |T(z)p = 

l-|i?+(z,t)|2 = l-|i?+(z,0)|2. 
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Moreover, we set 

Toizo) = T{0,zo) 



(2.9) 



= < 



f n i4r\ 

a-e(Co:0) 

n iarMexp(5i.7iog(|T(s)i)f y 
n iarMexpf5i./iog(ir(s)i)f 

.ae(-i,o)u(Co,i) / \ T 



2oe(-i,0), 

^oG(0,l), 



and 



az 



2 = 



f E (c^-a), 

ae(Co,o) 



(2.10) 



^oe(-i,o), 
kol = 1, 



E (c^-a) + ii/iog(|T(s)i)f, 

Cfce(-i,o) w 

E (C^-a) + ir/log(in5)|)fl, ^oe(0,i) 

Ufce(-i:0)u(Co,i) T 



In other words, T(z, zq) = To(zo)(l + Tx{zq)z + ©(z^)). 

Theorem 2.1 (Soliton region). Assume (|2.2p /or some ^ <E N and abbreviate by 

Ck = — 2io (ir II ^^e velocity of the k'th soliton determined by Re($(('/j)) = 0. Then 

the asymptotics in the soliton region, \n/t\ > 1 + C/ilog(t)^ for any C > 0, are as 

follows. 

Let e > sufficiently small such that the intervals [ck — £,Cfc + e], 1 < k < N , 
are disjoint and lie inside (— oo, —1) U (l,oo). 

// It ~ Cfcl < £ for some k, the solution is asymptotically given by a single soliton 



\{{2a{j,t))^To{zo 

oo ^ 

(2.11) Y. Kj,t)^-T,{ 



20 j 



l-Cfc+7fc(",^) 
l-Cfe+7fe(" + l,i) 

7fc(n,i)Cfe(l-Cfc) 



o{t- 



]=n+l 



2-^-' ' 2((7fe(n,i)-l)C,2 + l) 
where 

(2.12) 7fc(",i) = lkTiCk,-Ck - ^cl - l)-2e*K'=-';^")Cf 

// I J — Cfcl > £, /or aZ/ fc, one has 



o{t- 



(2.13) 



n(2a(j-,t))=To(zo)(l + 0(t-')) 

-n 
oo ^ 

^ 6(j,t) = -Ti(zo) + 0(i-')- 



j=n 



3=n+l 



Note that one can choose |j — Cfc| < ei for the regions where (|2.1ip is vahd, 
respectively | j — Cfc| > £2 for the regions where (|2.13p is vahd, such that the regions 
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overlap if ei > 62- Due to the exponential decay of the one-soliton solution, both 
formulas of course produce the same result on the overlap. 

In particular, we recover the well-known fact that the solution splits into a sum of 
independent solitons where the presence of the other solitons and the radiation part 
corresponding to the continuous spectrum manifests itself in phase shifts given by 
T(Cfc, —Ck — \/c-i — 1)~^- Indeed, notice that for Cfc G (—1, 0) this term just contains 
the product over the Blaschke factors corresponding to solitons Q with C,k < Q. 
For Cfc € (0,1) we have the product over the Blaschke factors corresponding to 
solitons Cj G (— 17 0), the integral over the full unit circle, plus the product over the 
Blaschke factors corresponding to solitons Q with (k > Cj ■ 

Furthermore, this result shows that in the region j > 1 the solution is asymp- 
totically given by a iV_-soliton solution, where N^ is the number of Q G (~ljO), 
formed from the data Q, 7j for all (k € ("l, 0). Similarly, in the region j < — 1 the 
solution is asymptotically given by a 7V-|_-soliton solution, where iV+ is the number 
of Cj G (0, 1), formed from the data Q, 7j for all Cj G (0, 1), where 

In the remaining region, we will show 

Theorem 2.2 (Similarity region). Assume (j2.2p with I > 5, then, away from the 
soliton region, \n/t\ < 1 — C for any C > 0, the asymptotics are given by 

1/2 



(2.15) +o(r")). 



j=ra+l ^ \ 0) / 

(2.16) +0(t-"), zo = e''°, 
for any a < 1. Here 

K^o)=--log(|T(zo)|), 
n 

$o(zo)-2(sin(0o)-^ocos(eo)), 

(2.17) Sizo) ^n/A - Si^izo) log |2sin(0o)| + 2arg(f(zo)) - arg(i?+(zo,0)) 

+ arg(r(it/(zo))), 

-1 / ^^ 

ra(\n^ ^~^'= \27ri7 V|r(2;o)|/s-zo s 

and r(z) is the gamma function. 

For a{n,t) respectively b{n,t) we obtain as a simple consequence: 
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Figure 3. Numerically computed solution a{n, 150) of the Toda 
lattice in Flaschka's variables. 

Corollary 2.3. Under the same assumptions as in Theorem \2.2\ we have 

a{n, t) =i + ('z£!!l(|Mf^') cos (i$o(zo) + ^{z^) \og{t) - <5(zo) - ^o) 

(2.18) +0(t-°), 

(2.19) +0(i""). 

Proof. To get the first formula for we use a{n, t) — \ n?ln(2'^0' ^))/ Iljln+i (^^(jj ^))- 
Now set x = f and observe O^m^) ^ 9q{x+\) ^ eQ{x)±e' {x)\+0{t-'^) uniformly 
in |a;| < 1 — C . Similarly for for the other terms and hence on checks that the only 
difference up to 0{t^°') errors in the above formulas for n and n ± 1 is a ^26*0 in 
the argument of the cosine (stemming from the t<I>o(zo) term). The second formula 
follows in the same manner from h{n, t) — X^i^n ^0' ^) ~ X]i=n+i ^ij^ 0- I— ' 

This is illustrated in Figure [Sj which shows the same solution as in Figure [2] but 
in Flaschka's variables. It is also interesting to look at the relation between the 
energy A of the underlying Lax operator H and the propagation speed at which the 
corresponding parts of the Toda lattice travel, that is, the analog of the classical 
dispersion relation. By the above theorems, the nonlinear dispersion relation is 
given by (see Figure |4|) 

(2.20) v{\) = ^, 
where 
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Figure 4. Nonlinear dispersion relation for the Toda lattice. 

We will not address the asymptotics in the missing region around \n\ « t. In the 
case \R+{z, 0)| < 1 the solution can be given in terms of Painleve II transcendents. 
If |i?+(z,0)| — 1 (which is the generic case), an additional region, the coUisionless 
shock region, will appear where the solution can be described in terms of elliptic 
functions. For the Painleve region we refer to [4], [20]. For the coUisionless shock 
region an outline using the ^-function method was given in |10j (for the case of the 
Korteweg-de Vires equation). The case of the Toda lattice will be dealt with in 

m- 

We also remark that the present methods can also be used to obtain further 
terms in the asymptotic expansion [S]. 

Finally, note that one can obtain the asymptotics for n > from the ones for 
n < by virtue of a simple reflection. Similarly for t > versus t < 0. 

Lemma 2.4. Suppose a{n,t), b(n,t) satisfy the Toda equation (j2.ip . then so do 

d{n, t) — a{—n — 1, i), &(n, t) — —b{—n, t) 

respectively 

a{n, t) — a{n, — t), 6(n, t) — —h{n, —t). 

3. The Inverse scattering transform and the Riemann-Hilbert 

PROBLEM 

In this section we want to derive the Riemann-Hilbert problem from scattering 
theory. The special case without eigenvalues was first given in Kamvissis [20] . How 
eigenvalues can be added was first shown in Deift, Kamvissis, Kriecherbauer, and 
Zhou [TTj. We essentially follow [5S] in this section. 

For the necessary results from scattering theory respectively the inverse scatter- 
ing transform for the Toda lattice we refer to [33] , [37] , [3H] . 

Associated with a{t),b{t) is a self-adjoint Jacobi operator 

(3.1) H{t)=a{t)S+ +a-{t)S- +b{t) 

in ^^(Z), where S^f{n) = f^{n) = /(n± 1) are the usual shift operators and £^(Z) 
denotes the Hilbert space of square summable (complex-valued) sequences over Z. 
By our assumption (|2.2p the spectrum of H consists of an absolutely continuous part 
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[—1,1] plus a finite number of eigenvalues Afe G M\[— 1, 1], 1 < fc < A^. In addition, 
there exist two Jost functions ip±{z,n,t) which solve the recurrence equation 

(3.2) iJ-(OV± {z, n,t)^ ^ '^^ V± [z. n, t), \z\<l, 
and asymptotically look like the free solutions 

(3.3) lim z=F'V±(z,n,t) = 1. 

n — ►ioo 

Both 'ip±{z,n,t) are analytic for < |z| < 1 with smooth boundary values for 
\z\ — 1. The asymptotics of the two Jost function are 

(3.4) i;±{z,n,t) ^ [I + 2B±{n,t)z + 0{z^)j, 
as z — > 0, where 

CXD CSO 

A+{n,t) = l[2a{j,t), B+{n,t) = - ^ b{j,t), 

j—n 7— n+1 

3.5 ' ' 

^ ^ n— 1 n— 1 

A-{n,t)= n MJ,t), B^{n,t)^- ^ 6(j,t). 

J — — 00 J — — 00 

One has the scattering relations 



(3.6) T{z)'ip^{z,n,t) = ip±{z,n,t) + R±{z,t)ip±{z,n,t), \z\ = I, 

where T{z), R±{z,t) are the transmission respectively reflection coefficients. The 
transmission and reflection coefficients have the following well-known properties 
([an Sect. 10.2]): 

Lemma 3.1. The transmission coefficient T{z) has a meromorphic extension to 
the interior of the unit circle with simple poles at the images of the eigenvalues Q ■ 
The residues ofT{z) are given by 

(3.7) Resc, T{z) = -Ck^^^^ = -Cfe7-,fe(i)Mfc(i), 

where 

(3.8) 7±Atr' ^T.\^±(^k,n,t)\^ 

and ip^{(k,n,t) = ^fe(OV'+(Cfc, ", 0- 
Moreover, 



(3.9) T{z)R+{z,t) + T{z)R^{z,t)^0, |T(z)|^ + |i?±(z, t)|^ = 1. 

In particular one reflection coefficient, say R{z,t) = R^{z,t), and one set of 
norming constants, say jk{t) = J+,k{t), suffices. Moreover, the time dependence is 
given by ([381 Thm. 13.4]): 

Lemma 3.2. The time evolutions of the quantities R-\^{z,t), 7+^fc(t) are given by 

(3.10) i?(z,i) = i?(z)e*("-""') 

(3.11) ^k{t)^lke"^^''-^^'\ 
where R{z) — R{z, 0) and jk — 7/c(0). 
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Now we define the sectionally meromorphic vector 

n 12) n^iz r, n - / {T{z)^l^^{z,n,t)z- i,+ {z,n,t)z-) , \z\ < 1, 

[6Al) ™^^'"'^^-| (Y,^(;2-i,„,i)2" r(z-i)^_(z-i,n,i)z-"), |z|>l. 

We are interested in the jump condition of m{z,n,t) on the unit circle T (ori- 
ented counterclockwise). To formulate our jump condition we use the following 
convention: When representing functions on T, the lower subscript denotes the 
non-tangential limit from different sides, 

(3.13) m±(z)- lini m(C), \z\ = 1. 

C^z.. |CI±i<i 

In general, for an oriented contour S, m^{z) (resp. m^{z)) will denote the limit 
of to(C) as C, —> z from the positive (resp. negative) side of S. Here the positive 
(resp. negative) side is the one which lies to the left (resp. right) as one traverses 
the contour in the direction of the orientation. Using the notation above implicitly 
assumes that these limits exist in the sense that m{z) extends to a continuous 
function on the boundary. 

Theorem 3.3 (Vector Riemann-Hilbert problem). Let S+{H{0)) = {R{z), \z\ = 
1; (Cfci7fc)j 1 < fc < N} the left scattering data of the operator H{0). Then m{z) — 
m{z, n, t) defined in p.l2p is a solution of the following vector Riemann-Hilbert 
problem. 

Find a function m{z) which is meromorphic away from the unit circle with simple 
poles at (k, Ck '^'^'^ satisfies: 
(i) The jump condition 

(3.14) 7n+(z) = ?7i_(z)u(z), i7(z) 



l-|i?(z)|2 -i?(z)e-**(^) 
i?(z)e**(^) 1 



for z G T, 
(ii) the pole conditions 



Res,-i m(z) = lim m(z) _ '^ ' 
'»'= z^c^^ \0 

(iii) the symmetry condition 

(3.16) miz-')^m{z)(j J 

(iv) and the normalization 

(3.17) m{0) = {mi 7712), tiii ■ m2 — 1 mi > 0. 

Here the phase is given by 

n 

(3.18) ^{z)^ z- z-'^ + 2-\ogz. 

Proof. The jump condition p.l4p is a simple calculation using the scattering re- 
lations (|3.6p plus (|3.9p . The pole conditions follow since T{z) is meromorphic in 
|z| < 1 with simple poles at Cfe and residues given by (j3.7p . The symmetry condition 
holds by construction and the normalization p.l7p is immediate from the following 
lemma. D 
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Observe that the pole condition at (k is sufficient since tlie one at (^^ follows 
by symmetry. Moreover, it can be shown that the solution of the above Riemann- 
Hilbert problem is unique '25'. However, we will not need this fact here and it will 
follow as a byproduct of our analysis at least for sufficiently large t. 

Moreover, we have the following asymptotic behaviour near z — 0: 

Lemma 3.4. The function m{z,n,t) defined in (I3.12p satisfies 

(3.19) m{z,n,t)=(^A{n,t){l-2B{n-l,t)z) -^^{l + 2B{n,t)z)^ + 0{z^). 

Here A{n,t) = A+(ri,i) and B(n,t) = Bj^{n,t) are defined in p.Sp . 

Proof. ThisfoUowsfromdSll) andr(z) ==yl+A_(l-2(B+-& + B_)z + 0(z2)). Q 

For our further analysis it will be convenient to rewrite the pole condition as a 
jump condition and hence turn our meromorphic Riemann-Hilbert problem into a 
holomorphic Riemann-Hilbert problem following llj. Choose e so small that the 
discs |z — Cfcl < e are inside the set {z|0 < \z\ < 1} and do not intersect. Then 
redefine ttj in a neighborhood of (k respectively C^ according to 



(3.20) m(z) 




\z-Ck\ < e, 
\z-'~Ck\<e, 



else. 

Then a straightforward calculation using Res,;; m = lim2^^(z — C)m(z) shows 

Lemma 3.5. Suppose m{z) is redefined as in p.20p . Then m{z) is holomorphic 
away from the unit circle and satisfies (|3.14p . p.l6p . p.l7p and the pole conditions 
are replaced by the jump conditions 

0\ 

m+{z) = m_(z) C)=7fce'*<«fc) -, h \z-(k\^e, 



(3.21) 




1 



TO+(z)=m_(z)( 2-C, 1, \z^-Ck\=s, 



where the small circle around (k *s oriented counterclockwise and the one around 
C^ is oriented clockwise. 

Finally, we note that the case of just one eigenvalue and zero reflection coefficient 
can be solved explicitly. 

Lemma 3.6 (One soliton solution). Suppose there is only one eigenvalue and a 
vanishing reflection coefficient, that is, S+{H{t)) = {R{z) = 0, \z\ ~ 1; (C,7)} with 
C e (-1,0) U (0,1) and 7 > 0. Then the Riemann-Hilbert problem PTTI)) - ([5T71) 
has a unique solution is given by 

(3.22) 

mo{z) = {f{z) f{l/z)) 



Vl - C + 7(^, Vl - C + C^lin, t) 



. z 
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where ^{n,t) — 76**'^'. In particular, 



(3.23) A+(n,t} = J -5- — - — -—2, B+{n,t} - 



i-e+jin,t)e' +^ ' ^ 2{i- e+^{n,t)ey 

Furthermore, the zero solution is the only solution of the corresponding vanishing 
problem where the normalization is replaced by m(0) — (0 7712) with TO2 arbitrary. 

Proof. By symmetry, the solution must be of the form mo{z) = [f{z) f{l/z)), 
where f{z) is meromorphic in C U {c)o} with the only possible pole at C. Hence 



A\ z-C, 

where the unknown constants A and B are uniquely determined by the pole con- 
dition ReSij /(z) — —C,^{n,t)f{(l^^) and the normalization /(0)/(oo) = 1, /(O) > 
0. D 

4. Conjugation and deformation 

This section demonstrates how to conjugate our Riemann-Hilbert problem and 
deform the jump contours, such that the jumps will be exponentially close to the 
identity away from the stationary phase points. In order to do this we will assume 
that R{z) has an analytic extension to a strip around the unit circle throughout 
this and the following section. This is for example the case if the decay in (j2.2|) is 
exponentially. We will eventually show how to remove this assumption in Section [S] 

For easy reference we note the following result which can be checked by a straight- 
forward calculation. 

Lemma 4.1 (Conjugation). Assume that ECS. Let D be a matrix of the form 

U{z)-^ 



(4-1) D{z)-^ Q ^(^, 

where d : C\S -^ C is a sectionally analytic function. Set 

(4.2) to(z) = m{z)D{z), 
then the jump matrix transforms according to 

(4.3) v{z) = D^{z)-\{z)D+{z). 

If d satisfies d{z^^) = d{z)^^ and d{0) > 0. Then the transformation rh{z) = 
m[z)D{z) respects our symmetry, that is, m{z) satisfies (j3.16p if and only if m{z) 



z e E\E, 



ze E. 

In order to remove the poles there are two cases to distinguish. If Re(<I>(CA;)) < 
the corresponding jumps p.2ip are exponentially close to the identity as t -^ 00 
and there is nothing to do. Otherwise, if Re(<I'(Cfe)) < 0, we use conjugation to 



does. 








In particular. 


we 


obtain 




(4.4) 




V- ( "" 
\v21d ^ 


Vl2d'\ 
W22 J ' 


respectively 








(4.5) 




/ d- 

V = 1 

\v2id^ d_ 


Vi2d+d. 

I rf+ 
3^ "22 
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turn the jumps into exponentially decaying ones, again following Deift, Kamvissis, 
Kriecherbauer, and Zhou [11] (see also [25|)- For this purpose we will use the next 
lemma which shows how 7fce*^''''^ can be replaced by its inverse. It turns out that 
we will have to handle the poles at C^; and (^^ in one step in order to preserve 
symmetry and in order to not add additional poles elsewhere. 

Lemma 4.2. Assume that the Riemann-Hilbert problem for m has jump conditions 
near ^ and i^^^ given by 



(4.6) 



m+(z) =TO_(z) ( _y^ 
m+{z) = TO_(z) ( ^ 



fZ \ 


k-Cl=£, 


-r) 


\z-'-C\=e 



Then this Riemann-Hilbert problem is equivalent to a Riemann-Hilbert problem for 
rh which has jump conditions near C and C,^^ given by 



m^{z) = TO_(z) 



1 



C(z-C)7 

1 



m+(z)=TO_(z)[ (^_t^)2 , I 

\Cz(Cz-i)t V 

and all remaining data conjugated (as in Lemma\4.1\l by 



\z-C\=e, 



(4.7) 



D{z) = ^^1 





Proof. To turn 7 into 7 ^, introduce D by 



D{z) 



and note that D{z) is analytic away from the two circles. Now set m(z) — 
m{z)D{z), which is again symmetric by D{z^^) = ( ? o)^(-z)( ? J)- The jumps 
along |z — CI = e and \z~^ — C,\ = e follow by a straightforward calculation and the 
remaining jumps follow from Lemma |4. II D 

The jumps along T are of oscillatory type and our aim is to apply a contour 
deformation which will move them into regions where the oscillatory terms will 
decay exponentially. Since the jump matrix v contains both exp(t$) and exp(— t<i>) 
we need to separate them in order to be able to move them to different regions of 
the complex plane. For this we will need the following factorizations of the jump 
condition ([531]) . First of all 



(4.8) 



v{z)^b-{z)-^b+{z), 



where 

6_ 
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This will be the proper factorization for z > zq. Here z > zq has to be understood 
as A(z) > Aq. Similarly, we have 

(4.9) viz) = B^iz)-' ( ' ^ ^' 1 p+(^), 



where 



1 0\ /i fl(^)c''^<'' 



V i-ifi(^)p V V" 1 / 

This will be the proper factorization for z < zq. 

To get rid of the diagonal part we need to solve the corresponding scalar Riemann- 
Hilbert problem. Moreover, for zo G (—1,0) we have Re($(2r)) > for z e (CoiO) 
and Re($(z)) < for z e (-l,Co) U (0, 1), for zq e T we have Re($(z)) > for 
z e (-1, 0) and Re($(z)) < for z e (0, 1), and for zq G (0, 1) we have Re($(z)) > 
for z e (-1,0) U (Co,l) and Re($(z)) < for z e (0,Co) (compare Figure [5] and 
note that by Re(<i>(z~^)) = — Rc($(z)) the curves Re($(z)) = are symmetric 
with respect to z i-^ z^^)- 

Together with the Blaschke factors needed to conjugate the jumps near the eigen- 
values, this is just the partial transmission coefficient T(z,zo) introduced in (|2.8p . 
In fact, it satisfies the following scalar meromorphic Riemann-Hilbert problem: 

Lemma 4.3. Set E(zo) = for zq G (-1, 0), I](zo) = {z e T| Re(z) < Re(zo)} for 
Zq G T, and S(zo) = T for Zq G (0,1). Then the partial transmission coefficient 
T{z, Zq) is meromorphic for z G C\I](zo), with simple poles at Qj and simple zeros 
at C~ for all j with i((^j + (^ ) < Aq, and satisfies the jump condition 

T+{z, zo) = r-(z,zo)(l - |i?(z)n, z G S(zo). 

Moreover, 

(i) r(z-i,z o)=r(z,zo)-\ zgC\S(zo), andr(0,zo)>0, 
(ii) T(z,zo) = T{'z,zq), z G C and, in particular, T{z,zq) is real-valued for 

zgR, 
(iii) r(z, zo) = r(z)(C + o(l)) with C ^ for \z\ < 1 near ±1 if ±1 G T,{zo) 
and continuous otherwise. 

Proof. That Q are simple poles and C~^ are simple zeros is obvious from the 
Blaschke factors and that T{z, zq) has the given jump follows from Plemelj's for- 
mulas. (i)-(iii) are straightforward to check. D 

Observe that for Co < Cn if Cn G (0, 1) respectively Co < 1 else we have T{z) = 

r(z,zo). 

Moreover, note that (i) and (ii) imply 
(4.10) |T(z,zo)P=r(z,zo)r(z,zo)=r(z-\zo)T(z,zo) = l, z G T\S(zo). 
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Now we are ready to perform our conjugation step. Introduce 



D{z) = { 



where 



<^-^ \d^{z), |2-i-a|<e, a, <i(Co+Co"'), 










Do{z) = 



T{z,zo)-' 




else, 





T{z,zo) 



Note that we have 



Now we conjugate our vector m(z) defined in p.l2p respectively (|3.20p using D{z), 

(4.11) m{z) = m{z)D{z). 

Since, by Lemma 14.31 (iii), r(z,zo) is either nonzero and continuous near z = ±1 
(if ±1 ^ 2](zo)) or it has the same behaviour as T{z) near z = ±1 (if ±1 G E(zo)), 
the new vector Tri{z) is again continuous near z = ±1 (even if T{z) vanishes there). 
Then using Lemma l4. II and Lemma 14.21 the jumps corresponding to eigenvalues 
\k <\{Co + Qo^) (if any) are given by 

f)(z)=r a7.T(.!.or^c-K.) J ^ \z-Ck\=e, 

\CkzikT(z,zoyc*'^«k) -^y 

and corresponding to eigenvalues Afe > ^(Co + CcT^) (if any) by 



(4.13) 



v{z) 



v{z) 



1 






1 



vO 



1 



k-Cfe| =e, 



a I 



In particular, an investigation of the sign of Re(<I'(z)) (see Figure [5] below) shows 
that all off-diagonal entries of these jump matrices, except for possibly one if Cfco — 
Co for some fcp, are exponentially decreasing. In the latter case we will keep the 
pole condition for ^fe,, = ^q which now reads 

(4.14) 



Res^-i rh(z) = lim rh(z) 



z^C-- 



0.o7^(4o,^o)-'e**(C.„) 




Furthermore, the jump along T is given by 

'6_(z)-i6+(z), A(z) > Ao, 
B_(z)-iB+(z), \{z) < Ao, 



(4.15) 



w(z) 
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^oe(-i,0) 



zoGT 




ZQ e (0, 1) 




Figure 5. Sign of Rc(<i>(z)) for different values of zq 



where 
(4.16) 
and 

(4.17) 



b-{z) = 



T(z 



1 



b+{z) = fl(z)e"^<'' 



T{z,zor 



B-{z) = 



B+iz) 



T.(z,zo)-\ 
l~R{z)R{z-^} 



R{z)e 



t<i>(z) 



.0 



T+{z,zo) 



i-R{z)R{z-^)^v-' ;•- 



-l\p-t$(2 



1 



Here we have used T±{z~^ , zq) = T±{'z, zq) = T±{z, zq) and R{z^^) = R{z) — R{z) 
for z S T to show that there exists an analytic continuation into a neighborhood of 
the unit circle. Moreover, using 

T±(z,zo) =T^{z''^,zo)~'^, zei;(zo), 

we can write 

.^ -^gx r-(z,zo)-^ ^ T_{z,zo) T+{z,zo)^ ^ T'+(z,zo) 

^ ■ ^ 1 - i?(z)i?(z-i) T_(z,zo)' 1 - i?(z)i?(z-i) T+(z,zo) 

for z g T, which shows that the matrix entries are in fact bounded. 

Now we deform the jump along T to move the oscillatory terms into regions 
where they are decaying. There are three cases to distinguish (see Figure [^1): 

Case 1: zq € (—1,0). In this case we will set E± = {z\ \z\ = (1 — e)*^} for some 
small e G (0, 1) such that E± lies in the region with ±Re($(z)) < and such that 
we do not intersect the original contours (i.e., we stay away from C )• Then we 
can split our jump by redefining m{z) according to 



(4.19) 



77l(z) 



(1-£)<|Z|<1, 

l<|z|<(l-e)-i, 
else. 



to(z)6+(z) 1, 
■m{z)b-{z)^^, 
m{z), 

It is straightforward to check that the jump along T disappears and the jump along 
E± is given by 

b+{z), z G S+, 

L(z)-i, zes_. 



(4.20) 



v{z) 



The other jumps (|4.12p . (|4.13p as well as the pole condition (|4.14p (if present) are 
unchanged. 
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Rc(*)<0 



Rc(*)>0 



Figure 6. Deformed contour 



Note that the resulting Rieniann-Hilbert problem still satisfies our symmetry 
condition ()3.16|) since we have 



By construction all jumps (14.121) . ()4.13p . and ()4.19|) are exponentially close to the 
identity as t — > oo. The only non-decaying part being the pole condition (I4.14p (if 
present). 

Case 2: zq G T \ {±1}. In this case we will set I]± — S]j_ U Ej. as indicated in 
FigurelH Again note that Ej_ respectively EH- lies in the region with ± Re($(z)) < 
and must be chosen such that we do not intersect any other parts of the contour. 
Then we can split our jump by redefining rh{z) according to 





rh(z)b^{z) ^, z between T and E?,_, 




to(z)6_(z)-\ z between T and EL, 


(4.21) m(z) = < 


m{z)B+{zy^, z between T and E^, 




rhiz)B^{z)-\ z between T and T,l , 




fn{z), else. 


One checks that the jump along T disappears and the jump along E± is given by 


(b+{z), zeEi., 


*<='^|^S"' ::S; 




[s_(z)-i, zeE2_. 



All other jumps (|4.12p and (|4.13p are unchanged. Again the resulting Riemann- 
Hilbert problem still satisfies our symmetry condition (|3.16p and the jump along 
E± away from the stationary phase points zq, Zq is exponentially close to the 
identity as i — > oo. 

Case 3: zq & (0, 1). In this case we will set E± = {z\ \z\ = (1 — e)"*"^} for some 
small e e (0, 1) such that E± lies in the region with q= Re($(z)) < and such that 
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we do not intersect the original contours. Then we can spht our jump by redefining 
m(z) according to 

(rhiz)B+{z)-\ (l-£)<|z|<l, 

(4.23) m{z)^ lm(z)B^{z)-\ 1< |z| < (1 - e)-i, 

[m(z), else. 

One checks that the jump along T disappears and the jump along S± is given by 

(4.24) v{z) = 

The other jumps (|4.12p . (|4.13p as well as the pole condition (|4.14p (if present) 
are unchanged. Again the resulting Riemann-Hilbert problem still satisfies our 
symmetry condition (|3.16p and all jumps (|4.12p . (|4.13p . and (|4.23p are exponentially 
close to the identity as t ^ cx3. The only non-decaying part being the pole condition 
(|4?T4l) (if present). 

In Case 1 and 3 we can immediately apply Theorem IB.6I to m as follows: If 
I7 — Cfel > £ for all fc we can choose 70 = 0. Since the error between w* and Wq 
is exponentially small, this proves the second part of Theorem 12.11 in the analytic 
case upon comparing 




(4.25) m{z)=m{z) , ^ ^,_ . ^ 



T{z, zo) 

T(z,zo)- 



with (j3.19p . The changes necessary for the general case will be given in Section [6l 
Otherwise, if | j — Cfc| < £ for some k, we choose 79 = Jkin, t). Again we conclude 
that the error between w* and Wq is exponentially small, proving the first part of 
Theorem 12.11 The changes necessary for the general case will also be given in 
Section [Sj 

In Case 2 the jump will not decay on the two small crosses containing the sta- 
tionary phase points zq and Zq . Hence we will need to continue the investigation 
of this problem in the next section. 

5. Reduction to a Riemann-Hilbert problem on a small cross 

In the previous section we have shown that for zq G T\{±1} we can reduce 
everything to a Riemann-Hilbert problem for m(z) such that the jumps are of 
order 0{t~^) except in a small neighborhoods of the stationary phase points zq 
and Zq^. Denote by S'-^(zq ) the parts of E_|_ U E_ inside a small neighborhood of 
Zq. In this section we will show that everything can reduced to solving the two 
problems in the two small crosses T,'-'{zo) respectively Y,^{zq^). 

It will be slightly more convenient to use the alternate normalization 

(5.1) m(z) = -.m(z), A ^ TqA, 

such that 

(5.2) m(0) = (l j,). 

Without loss of generality we can also assume that E consists of two straight lines 
in a sufficiently small neighborhood of zq. 
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We will need the solution of the corresponding 2x2 matrix 



(5-3) 



M^iz) = M^(z)u(z), z e E*^, 



Af'^(cx)) =1, 

where the jump v is the same as for m(z) but restricted to a neighborhood of one 
of the two crosses E*^ = (S+ U S_) n {z\ \z — zq\ < e/2} for some small e > 0. 
As a first step we make a change of coordinates 



/r A\ /■ V-2sin(6'o , . ZQi 

2^01 V-2sm(6'o) 

such that the phase reads <i>(z) — i$o + ^C^ + 0{C^). Here we have set 

respectively cos(0o) — ~n/t, which implies 

$0 - 2(sin(0o) - do cos(0o)), ^"(^o) - 2ie-2'^° sin(0o). 

The corresponding Riemann-Hilbert problem will be solved in Section [A) To ap- 
ply this result we need the behaviour of our jump matrices near zq, that is, the 
behaviour of T{z, zq) near z — > zq. 

Lemma 5.1. Let zq £ T, then 

(5.5) T(z,zo)=(-z^ zz ) r(z,zo) 

V z- zoj 

where v — —— log(|r(zo)|) and the branch cut of the logarithm used to define z"' = 
gijyiog(z) jg chosen along the negative real axis. Here 

rhi \ TT i/-i^-Cfc"' (^ 7, (\T{s)\\s + zds 

T{z,zo)= [[ a ^ • exp h— / log ^-f^ 

,^|t„s 2-C/c I27ri7 \\T{zo)\J s-z s 

is Holder continuous of any exponent less than 1 at z = zq and satisfies T{zo, zq) G 
T. 

Proof. This follows since 

The property T{zf), zq) e T follows after letting z ^ zq in (|4.10p . D 

Now if z(C) is defined as in (|5.4p and < a < 1, then there is an L > such 
that 

|T(z(C),zo)-rr(zo,zo)e-i-'°s(-2-"(«o))|<L|Cr, 

where the branch cut of C,^'^ is tangent to the negative real axis. Clearly we also 
have 

ii?(z(c))-i?(zo)i<iicr 

and thus the assumptions of Theorem lA.ll are satisfied with 
r = i?(zo)f (zo, 2o)-2e3-iog(-2si„(eo)) 
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and the solution of (j5.3|) is given by 

zq Mo 



(5.6) M'^(z) =1- 



O 



1 



(-2sin(6lo)i)i/2z-zo \t^ 

'°-[j3 j ' 
(5.7) p = y^e'('^/'^-^'s(i?(^o))+arg(r(i^))) ^_2 sin(6'o))"^'''r(zo, zo)2c-'**«r 



Mn 



where 1/2 < a < 1, and cos(0o) = — Aq. Note |r| = |i?(zo)| and hence u = 
-J-logil-\R{zo)n. 

Now we are ready to show 

Theorem 5.2. The solution 7fi(z) is given by 

1 



(5.8) m{z) = (1 l) 



i-2sm{eoW/ 



YHimoiz) +mo{z)) + 



whe 



(5.9) mo{z)^[Pji^^ -/3l^)> mo{z) = moiJ)=mo{z-')(l j) 

Proof. Introduce m{z) by 



m{z) — < 



'm{z)M^{z)-\ \z-zo\<e, 
7h{z)M^{z)-'^, \z~'^-zo\<e, 
rh(z), else, 



where 



-i^^O 1 



Mn 



oil 



^' V Oj ^ V ^J {-2sm{9Q)ty/^ z- zq \t 

The Riemann-Hilbert problem for m has jumps given by 



{M^{z)-\ \z-zo\=e, 

M'^{z)v{z)M'^{z)'-\ ze±,§ <\z-zo\<s, 

I, ze S, |z-zo| < f, 

v{z) ^ i M^ {z)-\ |z-i-zo|=e, 

M^{z)v{z)M^{z)~\ z e E, I < |z-i - zol < £, 

I, ze Sjz^i-zol < |, 

,-0(0), else. 



The jumps are I + 0{t ^/^) on the loops 



2; - Zol = e, |z" 



Zol = e and even 



I + 0{t ") on the rest (in the L"^ norm, hence also in the L^ one). In particular, 
as in Lemma lA. 31 we infer 



||a*~(1 l)h = 0{t-'/^). 



24 H. KRUGER AND G. TESCHL 



Thus we have with Jloo as in (jB.SP 

to(z) = (l l) + 77^ / IJ.{s)w{s)nao{s,z) 



27ri js 
+ T^ / m(s)(M^(.s)-i _ l)n^(s, z) + o{t-^) 

(1 1) + , , w^. uu/2 (1 ^)Mo^f ^^n^^is,z) 

{~2sm{9o)ty/^ "■ 27ri 7|^_3„|^^ s - zo 

1 . x^r^ 1 /■ S 



(-2sin(6'o)t)i/^ ' 27ri7|^-i_3^l^^ s-zo 

= (^ ^)- (-2sin(W/- ("^°(^) + '^°(^» + ^(F) 
finishing the proof. D 

Hence, using ([3?T9| and ((5TT|l . 

(5.10) {rh{z))2 = 4 (1 + (^1 + 25)-~ + 0{z')) 

and comparing with 
(5.11) 

/ 2Re(/3) ^ ^ 2Re(z^/3) ^ ^ nr 2, ^ ^ H 



(-2 sin(0o)t)i/2 ; \^ (_2 sin(0o)i)i/ 

we obtain 

(512) .12 ,, 2Re(/3) /l 

^^•^^^ ^ "^+(-2sin(0o)t)V2+^l,i« 

and 

In summary we have 



which proves Theorem 12.21 in the analytic case. 

Remark 5.3. Note that, in contradistinction to Theorem \B.(A Theorem \5.^ does 
not require uniform boundedness of the associated integral operators, but only some 
knowledge of the solution of the Riemann-Hilbert problem. However, it requires that 
the solution is of the form I + o(l) and hence cannot be used in the soliton region. 
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6. Analytic Approximation 

In this section we want to present the necessary changes in the case where the 
reflection coefRcient does not have an analytic extension. The idea is to use an 
analytic approximation and split the reflection in an analytic part plus a small rest. 
The analytic part will be moved to the complex plane while the rest remains on 
the unit circle. This needs to be done in such a way that the rest is of 0{t^^) and 
the growth of the analytic part can be controlled by the decay of the phase. 

In the soliton region a straightforward splitting based on the Fourier series 

oo 

(6.1) R(z) = Y^ R(k)z^ 

k — — OG 

will be sufficient. It is well-known that our assumption (|2.2p implies k'-R{—k) e 
£i(N) (this follows from the estimate [381 eq. (10.83)]) and R e C'(T). 

Lemma 6.1. Suppose R{k) e t{I.), k^R{-k) e i^{N) and let < e < I, P > 
be given. Then we can split the reflection coejjicient according to R{z) = Ra,t{z) + 
Rr^t{z) such that Ra.t{z) is analytic in < \z\ < 1 and 

(6.2) \RaAz)e'^'\^0{t-'), l-£<|z|<l, \RrAz)\ ^ Oit-'), \z\ = 1. 

Proof. We choose Ra,t{z) — YlV=-K{t) R{k)z'^ with K{t) = L -io'^n~e') ^-l ^°^ some 
positive /3o < P. Then, for 1 — e < \z\, 



|i?a,t(z)e-^*| < y^ \R{k)\e-^\l^ef < ||i?||ie-^*(l-e)-^^(*) < \\R\\ie-'-^->^°^' 
f 

Similarly, for \z\ 



k=~K(t) 



^^*^-^ ^ fc' ,.. ,., const const 



\RrAz)\< l^ \R{k)\<con.st l^ |i?(_fc)|< < 

fc=-oo k=K(t) + l ^ ' ^ ' 



D 



To apply this lemma in the soliton region zg G ("IjO) we choose 
(6.3) /?= min -Re($(z))>0 

and split R{z) = Ra,t{z) + Rr.t{z) according to Lemma [6. II to obtain 

b±{z) = ba,t,±{z)hr,t,±{z) = br,t,±{z)ba.t,±{z). 

Here ba,t,±{z), br,t,±{z) denote the matrices obtained from b±{z) as defined in (|4.16p 
by replacing R{z) with Ra^z), Rr^tiz), respectively. Now we can move the analytic 
parts into the complex plane as in Section [4] while leaving the rest on T. Hence, 
rather then (j4.20p . the jump now reads 

{ba,t,+ {z), zeT,+ , 

ba,t,-{z)-\ zeE_, 

br,t,~iz)~^br,t, + iz), zeT. 

By construction we have v{z) = l+0(t^') on the whole contour and the rest follows 
as in Section m 
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In the other sohton region zq G (0, 1) we proceed similarly, with the only dif- 
ference that the jump matrices B±{z) have at first sight more complicated off 
diagonal entries. To remedy this we will rewrite them in terms of left rather then 
right scattering data. For this pm'pose let us use the notation Rr{z) = R+{z) 
for the right and Ri{z) = R^(z) for the left reflection coefficient. Moreover, let 
Tr{z,zo) = T{z,zo) be the right and Ti{z,zo) = T{z)/T{z, zq) be the left partial 
transmission coefficient. 

With this notation we have 



(6.5) v{z) 

where 



6_(z)-i6+(z), A(z) > Ao, 

B^{z)-^B+{z), A(z)<Ao, 






1 



r,(z,zo)2 



and 




\ \T{z] 

B+{z) = ^ \T{zW -"'■^^ '"^ 

Using (|3.9p together with (|4.18p we can further write 




Now we can proceed as before with B±{z) as with b±{z) by splitting Ri{z) rather 
than Rr{z). 

In the similarity region we need to take the small vicinities of the stationary 
phase points into account. Since the phase is quadratic near these points, we 
cannot use it to dominate the exponential growth of the analytic part away from 
the unit circle. Hence we will take the phase as a new variable and use the Fourier 
transform with respect to this new variable. Since this change of coordinates is 
singular near the stationary phase points, there is a price we have to pay, namely, 
requiring additional smoothness for R{z). We begin with 

Lemma 6.2. Suppose R{z) G C^iT). Then we can split R{z) according to 

(6.6) R{z)^Ro{z) + {z-zo){z-z^)H{z), z G S(zo), 

where Rq{z) is a real polynomial in z such that H[z) vanishes at zq^'Sq of order 
three and has a Fourier series 

oo 

(6.7) H{z)^ V iJ,e'="°*(^), u;o = J, , ^ , 

, 7rcos6io +*o 

with k Hk summable. Here ^o — ^(^o)/i- 
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Figure 7. Contours of a cross 

Proof. By choosing a polynomial Rq we can match the values of R and its first four 
derivatives at zqjJo. Hence H{z) e C"'(T) and vanishes together with its first three 
derivatives at zq^Jo. 

When restricted to T,{zo) the phase ^{z)/i gives a one to one coordinate trans- 
form S(zo) -^ [^^0, i*&o+ii^o] and we can hence express H(z) in this new coordinate. 
The coordinate transform locally looks like a square root near zo and 'zq, however, 
due to our assumption that H vanishes there, H is still C^ in this new coordinate 
and the Fourier transform with respect to this new coordinates exists and has the 
required properties. D 

Moreover, as in Lemma |6. II we obtain: 

Lemma 6.3. Let H{z) be as in the previous lemma. Then we can split H{z) 
according to H{z) — Ha.t{z) + Hr,t{z) such that Ha^t{z) is analytic in the region 
Re($(z)) <{) and 

(6.8) |7J,,,(z)e*(-)*/2| ^ 0(1), Re($(z)) < 0, \z\ < 1, \HrAz)\ = O(t-i), \z\ = 1. 

Proof. We choose Haj{z) = J2T=-K(t) Hke'"^'^^'"'' with K{t) = Li/(2^)J- The rest 
follows as in Lemma [01 D 



By construction i?a,i (2:) = Roiz) + {z— ZQ){z—Ja)Ha,tiz) will satisfy the required 
Lipschitz estimate in a vicinity of the stationary phase points (uniformly in t) and 
all jumps will be I + 0{t~^). Hence we can proceed as in Section[5l 

Appendix A. The solution on a small cross 

Introduce the cross S = Ei U • • • U S4 (see Figure [7]) by 

El = {uc"'''/^ u e [0, (»)} E2 = {ue''^/^ u e [0, 00)} 

(A.l) E3 = {ue3''^/^ u e [0, 00)} E4 = {^ie-3'^/^ u £ [0, 00)}. 

Orient E such that the real part of z increases in the positive direction. Denote 
by D = {z, \z\ < 1} the open unit disc. Throughout this section z"^ will denote 
the function e"''°s(^\ where the branch cut of the logarithm is chosen along the 
negative real axis (— (X),0). 
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Introduce the following jump matrices {vj for z e Sj) 

fl -i?i(z)z2'''e-**(^)\ f 1 

""'-[O 1 j ' ''^ - l^i?2(z)z-2-e**(^) 1 

Now consider the RHP given by 

(A. 3) m^{z) = m^{z)vj{z), z G S^, j = l,2,3,4, 

to(z) — > I, z — > CX). 

We have the next theorem, in which we follow the computations of Sections 3 and 
4 in [4]. The method can be found in earlier literature, see for example [19] , One 
can also find arguments like this in Section 5 in [20] or (3.65) to (3.76) in [6j. 
We will allow some variation, in all parameters as indicated in the next result. 

Theorem A.l. There is some po > such that Vj{z) — I for \z\ > po- Moreover, 
suppose that within \z\ < po the following estimates hold: 

(i) The phase satisfies $(0) = i$o G iK, *'(0) = 0, $"(0) = i and 

iz^ 
(A.5) |$(,)_$(0)_— |<C|z|^ 

(ii) There is some r G D and constants (a, L) € (0, 1] x (0,oo) such that Rj, 
i = 1, . . . , 4, satisfy Holder conditions of the form 

|i?i(2)-r|<L|zr, |i?2(z)-r|<L|zr, 

(A.6) |i?3(;^)„_!_|<L|zr, |i?4(^)__!_|<L|zr. 



Then the solution of the RHP (|A.3P satisfii 



es 



l + Q 



(A.7) ^(,) = i+i_i_(^o ^^j+or-^i 

for \z\ > Po, where 

(A.8) P = ^gi(^/4-arg(r)+arg(r(i^)))g-it$oi-ii'^ ^ = _ 1. log(l _ |^|2)^ 

27r 
Furthermore, if Rj{z) and $(z) depend on some parameter, the error term is uni- 
form with respect to this parameter as long as r remains within a compact subset 
of D and the constants in the above estimates can be chosen independent of the 
parameters. 

We remark that the solution of the RHP (|A.3P is unique. This follows from the 
usual Liouville argument [3l Lem. 7.18] since det(wj) = 1. 

Note that the actual value of po is of no importance. In fact, if we choose 
< pi < Po, then the solution rh of the problem with jump v, where v is equal to 
V for |z| < pi and I otherwise, differs from m only by an exponentially small error. 

This already indicates, that we should be able to replace Rj (z) by their respective 
values at z = 0. To see this we start by rewriting our RHP as a singular integral 
equation. We will use the theory developed in Appendix |B] for the case of 2 x 2 
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matrix valued functions with mo(z) = I and the usual Cauchy kernel (since we 
won't require symmetry in this section) 

ds 



S — z 

Moreover, since our contour is unbounded, we will assume w G L^C^) H i^(S). All 
results from Appendix |B] still hold in this case with some straightforward modifica- 
tions if one observes that /i — I e L^(I]). Indeed, as in Theorem IB. 3[ in the special 
case b+{z) ~ Vj{z) and b-{z) = I for z G Ej, we obtain 

(A.9) m{z) = I + — / ^,{s)w{s)^, 

27ri J^ s~ z 

where /Lt — I is the solution of the singular integral equation 

(A.IO) (I-C^)(/i-I) = C,,I, 

that is, 

(A.ll) M = I+(I-C^)"'Cu,I, C^f = C^{tvf). 

Here C denotes the usual Cauchy operator and we set w(z) = w+(z) (since w-^{z) = 

0). 

As our first step we will get rid of some constants and rescale the entire problem 
by setting 

(A.12) m{z) = D{t)-^m{zt-^''^)D{t), 

where 

(A.13) D{t)^(^^'l' ^0^), d(i)=e"*o/2i'^/^ d{t)-'=W)- 

Then one easily checks that rh{z) solves the RHP 

(A. 14) m+(z) = m_(z)%(z), z G E^, j = 1,2,3,4, 

m(z) — > I, z — »■ cxD, z ^Yi^ 

where %(z) = D{ty^Vj{zt''^l'^)D{t), j = 1, . . . ,4, explicitly 

'l -i?i(zt-l/2)^2i.g-t(*(^t-V2)_$(0))N 



"^1 (2) = , 

^'^^^ " U2(^i-'/')2-''V(*(-*"''')-*(")) ij ' 

„ ^ /l -i?3(zri/2)z2i.e-*(*(^*"'^')-*(")) 
^^^^'^ ^0 1 

(A.15) {;4(z) = (^^^(^^-i/2)^-2iv(*(-*-^^^)-*(0)) 1 

Our next aim is to show that the solution to(z) of the rescaled problem is close to 
the solution m^(z) of the RHP 

(A.16) rh%{z)^rff_{z)^{z), z € E^, j = 1,2,3,4, 

rn'^{z) ^ I, z ^ 00, z ^ E, 
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associated with the following jump matrices 

(A.17) «3W=(^o ^^, j, «4W = (^_^,-2i.,i.V2 ij- 

The difference between these jump matrices can be estimated as follows. 

Lemma A. 2. The matrices w'^ and w are close in the sense that 

(A.18) Wj(z) = u;^^(z) + 0(t-"/2e-l"l'/**), z € E^, j = l,...4. 

Furthermore, the error term is uniform with respect to parameters as stated in 
Theorem \A.1[ 

Proof. We only give the proof z e Si, the other cases being similar. There is only 
one nonzero matrix entry in Wj{z) — w'^{z) given by 

_ f-i?i(zi-i/2)z2-e-*(*(-*"''')-*(0)) +rz2-e-'^'/2, \z\ < p^t^/^, 

A straightforward estimate for \z\ < po^^ shows 
\W\ = e'''^/4|^^(^^-i/2)g-tl.(.t-V2) __|^„|,|2/2 

< e'''^/4|^^(^^-i/2) _ ^|e-NlV4 + e''^/H\Hzt-^/^)\e-\'\"'\ 

where 6(z) = $(z) - $(0) - ^z^ = ^—^z^ + Here we have used ^z^ = i|zp 

for z € El and Re(-i$(zi~i/2)) < \z\'^/4^ by jXI)) . Furthermore, by (fOj) and 

dm, 

IW'I < e'-^/^Li-^/^i^re-l^l'/^ + e''^/4Ct-i/2|^|3e-NlV4^ 

for |z| < pot^/'^. For |z| > po*^^^ we have 

1^1 < e^'^/^g-l^lVa < e'''^/4g-Pot/4g-|z|V4 

which finishes the proof. D 

The next lemma allows us, to replace m{z) by m'^{z). 
Lemma A. 3. Consider the RHP 

(A. 19) m+iz) = m^{z)v{z), z e E, 

m,{z) — > I, z ^ oo, z ^ E. 

Assume that w G ^^(E) n L°°(E). Then 

(A.20) ii/^-^ii^^r^^ 

provided c\\w\\oo < 1, where c is the norm of the Cauchy operator on ^^(E). 
Proof. This follows since /i = p — I G i^(E) satisfies (I — Cu])fi = C^I. □ 
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Lemma A. 4. The solution 771(2) has a convergent asymptotic expansion 

(A.21) m(z) = I+-M(i)+0(^) 

for \z\> pot^''^ with the error term uniformly in t. Moreover, 

(A.22) M{t) = Af= + 0(t-"/2). 

Proof. Consider rh'^{z) = m(z)m'^(z)^^, whose jump matrix is given by 

v'^iz) = m'L{z)v{z)v''{z)-^rh1{z)-^ = I + m'L{z){w{z) - w''{z))m''_{z)-\ 

By Lemma lA.21 we have that w — w'^ is decaying of order i^"' ^ in the norms of L^ 
and L°° and thus the same is true for w'^ — v'^ — I. Hence by the previous lemma 

Furthermore, by fi'^ = ml = m.^{rh'L)^^ = fi(Jl'^)^^ we infer 

||A-A1l2 = or"/2) 

since fi'^ is bounded. Now 

■m{z) — I — —^ — / fl{s)w{s)ds + -—^— / sfi{s)w{s) 

27ri z J-^ 27ri z J-^ s — z 

shows (recaU that w is supported inside \z\ < p^t^/'^) 

m(z)=I+iAf(t) + 0(Mf»!^<f^), 
z z^ 

where 

M{t) = ~^J fi{s)w{s)ds. 

Now the rest follows from 

M{t) = Ar - — / {fi{s)w{s) ~ p%s)w'{s))ds 

using Wfiw - fi^w^Wi < \\w -w^Wi + \\fi-l\\2\\w-w''\\2 + ||/i - Al|2||w'=||2. □ 

Finally, it remains to solve (jA.16[) and to show: 
Theorem A. 5. The solution of the RHP (|A.16[) is of the form 

(A.23) m^iz)=I+-AI' + Oi\), 

z z'^ 

where 

(A.24) ^^' = i(§ "0^)' /3 = y^e'('^/4-"sW+arg(rM))_ 

The error term is uniform with respect to r in compact subsets of J}. Moreover, the 
solution is bounded (again uniformly with respect to r). 

Given this result, Theorem I A. II follows from Lemma [A. 41 

m{z) = D{t)m{zt^'^)D{t)-^ = I + -^D{t)M{t)D{t)-^ + 0{z-h-^) 
(A.25) = I + -j^D{t)M^D{t)-^ + 0(f-(i+")/2) 

for l^l > po, since D{t) is bounded. 
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Figure 8. Deforming back the cross 



The proof of this result will be given in the remainder of this section. In order 
to solve ()A.16P we begin with a deformation which moves the jump to R as follows. 
Denote the region enclosed by M and Ej as ilj (cf. Figure [S]) and define 



(A.26) 
where 



Do{z), else, 



m (z) = TO (z) ' ^ ' ■> 



Do{z) 



z-'^-^e'"/^ I ' 



and 



Di = 



1 r 
1 



D2 = 



1 
r 1 



D.^ 



l-|rP 
1 



D4 = 



M^ 



Lemma A. 6. The function rh^{z) defined in (JA.26P satisfies the RHP 



(A.27) m^(z)=TOl(z) 



1 



l-lrl- 



TO'=(z) = (I+-Ar + ...)L>o(^), 



z e 



z^oo, - <arg(z) < — . 



Proof. First, one checks that to^(z) — ■m'L{z)DQ{z)^^vl{z)DQ(z)Di = rh'L{z), z £ 
El and similarly for z G E2,5]3,S4. To compute the jump along R observe that, 
by our choice of branch cut for z"', Dq{z) has a jump along the negative real axis 
given by 



DoA^) = 



,{log\z\±iTT)iiy^-iz^/4 







-{log\z\±i7T)iu iz^/4 



z < 0. 



Hence the jump along R is given by 

D^^D2, z>0 and D^^D^l{z)Do.,+{z)D3, z < 0, 



and (JA.27P follows after recalling e"^''" = 1 - 



D 
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Now, we can follow (4.17) to (4.51) in [4 to construct an approximate solution. 

The idea is as follows, since the jump matrix for (|A.27[) . the derivative ■j^rrf{z) 
has the same jump and hence is given by n{z)TTf'{z), where the entire matrix n{z) 
can be determined from the behaviour z — > cxd. Since this will just serve as a 
motivation for our ansatz, we will not worry about justifying any steps. 

For z in the sector ^ < arg(z) < ^ (enclosed by S2 and E3) we have 'nf{z) = 
m'^{z)Do{z) and hence 



( — to'=(z) + — cr3m'=(z) j m'=(z) ^ 



i^ z, ^ .., . d ^ ,., , .z ^ .. , . \ ^ ,. , , _i 



i(- - -)m=(z)cr3 + -ii^m^iz) + \-(j3,rrf{z) ) to'=(z) 

=^[.3,Mi+o(i), .3=(; \ 

Since the left hand side has no jump, it is entire and hence by Liouville's theorem 
a constant given by the right hand side. In other words, 

(A.28) A^^(z) + |a377V(z) ^ /3m^(z), /? = (^^^ ^^^^'-[a,,An 

This differential equation can be solved in terms of parabolic cylinder function 
which then gives the solution of (|A.27p . 

Lemma A. 7. The RHP (|A.27p has a unique solution, and the term, M'^ is given 
by 

(A.29) M- = i(^ ~'1'A, /3i2=^=V^e'('^/4-^''gW+^^g(r("'»). 

VP21 U J 

Proof. Uniqueness follows by the standard Liouville argument since the determinant 
of the jump matrix is equal to 1. To find the solution we use the ansatz 

yp2l{z) i^22{z)J 

where the functions ^Jjk{z) satisfy 

V'l'l(2) = - f 2 + 4^^ ~ ^12/321 j V'll(z), ^'12(2;) = -J— ( ^- YJ ^22{z), 
^2l{z) = ±-(^± + 11^ ^i,(z), V'^',(Z) = (^^ - iz^ + p^^p^^^ ^^^(^). 

That is, i^uie'^'^^^'^O satisfies the parabolic cylinder equation 

1 1 



with a — i/3i2/?2i and V'22(e"^/'*C) satisfies the parabolic cylinder equation with 
a = -i/?i2/32i- 
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Let Da be the entire parabolic cylinder function of §16.5 in [42, and set 

'/'iii^) - |e"'^/4A.(e'"/*^), Ini(z) < 0, 

f e'^'^/4i:)_i^(-ie'^/4z), Im(z) > 0, 

Using the asymptotic behavior 

i?,(z) = z«e-^'/4(l-^^|^ + 0(z-4)), z^oo, |arg(z)| < 37r/4, 
shows that the choice /3i2/32i — v ensures the correct asymptotics 
7^n(^) = z'-'e--'/^^^ ^ 0{z-^)), Vi2(^) = -i/?i2z-'V-'/*(z-^ + 0(^-')), 

^21 (^) = i/32iz'''e--'/4(z-i + 0(z-^)), fe(z) = z--e-'/4(l + 0{z-^)), 
OS, z -^ CO inside the half plane Ini(z) > 0. In particular, 



ra 



\z)^{^^-_^M^ + 0{z-^))Dq{z) with M^ = \i^_^^ M 



It remains to check that we have the correct jump. Since by construction both 
limits to'^(z) and ml (2;) satisfy the same differential equation (jA.28p . there is a 
constant matrix v such that ?7i^(z) — rh'L{z)v. Moreover, since the coefficient 
matrix of the linear differential equation (|A.28|) has trace 0, the determinant of 
TOj_(z) is constant and hence det(TOj_(2)) = 1 by our asymptotics. Moreover, a 
straightforward calculation shows 

V = ml(0)-im;(0) = ,^,.,/. -.„/. ^^^r^ 



v^r(-iiy) ' 



v^ _ /3: 



where 7 = -^ = ^. Here we have used 

' P12 ^Jv 



90/2 rz- 9(l+a)/2 rz 

plus the duplication formula r(z)r(z+i) = 2^^'^^ ^V{2z^ for the Gamma function. 
Hence, if we choose 

^ V2^e''^/4e-'^'^/2 *"' 

we have 

'\ - |r|2 -f 
r 1 

since hP = 1. To see this use |r(-iz/)P = r(i-it^)rM ^ . ?; which follows 

from Euler's reflection formula r(l — z)r(z) = ^^^7^^) for the Gamma function. 
In particular, 

which finishes the proof. D 
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Remark A. 8. An inspection of the proof shows that rh'^ is given by the solution of 
a differential equation depending analytically on v . Hence, rh'^ depends analytically 
on V = —-^XogiX ~ |rp). This implies local Lipschitz dependence on r as long as 

r e D. 

Appendix B. Singular integral equations 

In this section we show how to transform a nieromorphic vector Riemann-Hilbert 
problem with simple poles at C, C~^j 

m^{z) — m^{z)v{z), z G Ti, 

miz 1) =m(z) r ^j , 

m(0) = (1 ma) , 

where C G (—1,0) U (0,1) and 7 > 0, into a singular integral equation. Since 
we require the symmetry condition for our Riemann-Hilbert problems we need to 
adapt the usual Cauchy kernel to preserve this symmetry. Moreover, we keep the 
single soliton as an inhomogeneous term which will play the role of the leading 
asymptotics in our applications. 

Hypothesis B.l. Suppose the jump data (S,ti) satisfy the following assumptions: 

(i) S consist of a finite number of smooth oriented finite curves in C which 
intersect at most finitely many times with all intersections being transver- 
sal. 
(ii) E does not contain 0, C^^ . 

(iii) E is invariant under z i~> z^^ and is oriented such that under the mapping 
z i-^ z^^ sequences converging from the positive sided to E are mapped to 
sequences converging to the negative side. 

(iv) The jump matrix v is invertible and can be factorized according to v — 
bZ 6_|_ = (I — w_)^^(I + w+), where w± ~ ±(6-|- — I) are continuous and 
satisfy 



w±{z-')=(l J)^tW(? J) 



(B.2) w^{z-')=r Aw^iz)r M, zeS 



The classical Cauchy-transform of a function / : E — > C which is square inte- 
grable is the analytic function Cf : C\E ^ C given by 

(B.3) (^Cf){z) = ^ f ^^ds, zeC\E. 

Denote the non-tangential boundary values from both sides (taken possibly in the 
L^-sense — see e.g. [21 eq. (7.2)]) by C+/ respectively C_/. Then it is well-known 
that C+ and C_ are bounded operators -Zj^(E) -^ L'^(E), which satisfy C+ — C_ = I 
and C+(7_ = (see e.g. [Ij). Moreover, one has the Plemelj-Sokhotsky formula 

([S]) 

C± = i(ii/±I), 
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where 

(B.4) ^Hf){t)^--f i^ds, teS, 

is the Hilbert transform and ^ denotes the principal value integral. 

In order to respect the symmetry condition we will restrict our attention to the 
set is(S) of square integrable functions / : E — > C^ such that 



(B.5) 



/(=-)= /w(? ;) 



Clearly this will only be possible if we require our jump data to be symmetric as 
well (i.e., Hypothesis IB. II holds) . 

Next we introduce the Cauchy operator 

(B.6) iCf){z)^^ f f{.s)ncis,z) 

acting on vector- valued functions / : S ^ C^. Here the Cauchy kernel is given by 

(B.7) n^is,z)=(^^^ ^A'^'^i^'o^^ ^l^)d^' 

\ ^ s-Cs-zJ ^ s-z s-QJ 

for some fixed C, ^Yj. In the case C = oo we set 

(B.8) n^{s,z) = h-^-^ - _j_ 

\ s—z 

and one easily checks the symmetry property: 

(B.9) fic(lMlA)-(j J)f^c(^:^)(j I 

The properties of C are summarized in the next lemma. 

Lemma B.2. Assume Hvvothesis \B.l[ The Cauchy operator C has the properties, 
that the boundary values C± are bounded operators -^^(S) —)■ ^^(S) which satisfy 

(B.IO) C+-C_=I 

and 

(B.U) (C/)(ri) = (0 *), (C/)(C) = (* 0). 

Here * is a placeholder for an unspecified value. Furthermore, C restricts to Lg(E), 
that is 

(B.12) {Cf)iz-')^{Cf)iz)(j J), zeC\E 

for f e is(S) and if w± satisfy (|B.2[) we also have 

(B.13) C±{fw^){l/z)^C^{fw±)iz)(^^ J), zeE. 

Proof. Everything follows from (jB.9p and the fact that C inherits all properties 
from the classical Cauchy operator. D 
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We have thus obtained a Cauchy transform with the required properties. Fol- 
fowing Section 7 and 8 of [1] respectively j25], we can solve our Riemann-Hilbert 
problem using this Cauchy operator. 

Introduce the operator Cw ■ -^s(^) ^ ^s(^) by 

(B.14) C^f^C+{fw-)+C-{fw+), feLliY.) 

and recall from Lemnia l3.6l that the unique solution corresponding to u = I is given 

by 

mo(z) = (/(.) /(i)), f{z)-^^^^{ie'-^ + l-e) 

Observe that for 7 = we have /(z) — 1 and for 7 = 00 we have f{z) ~ Q"^ ^~f , . 
In particular, m{){z) is uniformly bounded away from C, for all 7 G [0, 00]. 
Then we have the next result. 

Theorem B.3. Assume Hvvothesis \B.l[ 

Suppose m solves the Riemann-Hilbert problem (IB. II) . Then 

(B.15) m(z) = (1 -co)mo(z) + — ^ / fj.{s){w+{s) + w-{s))Qi;{s, z), 



27ri ./s 



wh' 



ere 



/i = TO+fe I = rn^b_ and Cq ~ / /i(s)(w+(s) + ui_(s))51r(s, 0) 

V27riJE /i 

Here (rn)j denotes the j 'th component of a vector. Furthermore, /i solves 
(B.16) (I-C^)M=(l-co)mo. 

Conversely, suppose jl solves 
(B.17) (I-C^)A = mo, 

and 

co= (J-l^fi{s)i^+{s)+W-{s))nci'S,0)^ ^-1, 

then m defined via (jB.lSp . with (1 — cq) = (1 + cq)^^ anrf /i = (1 + co)^^/i, solves 
the Riemann-Hilbert problem (jB.ip anrf fi — m± b^ . 

Proof. If m solves (jB.l|) and we set fi ~ m±b^ , then m satisfies an additive jump 
given by 

771+ — r7i_ — /i(w+ + W-). 

Hence, if we denote the left hand side of (jB.15|l by rh, both functions satisfy the 
same additive jump. Furthermore, Hypothesis [BTT] implies that /i is symmetric and 
hence so is rh. Using (|B.11|1 we also see that rfi satisfies the same pole conditions 
as mo. In summary, m — m has no jump and solves (jB.ljl with 1; = I except for the 
normalization which is given by ?n(0) — m(0) = (O *). Hence Lemma 13.61 implies 
m — rh = 0. 

Moreover, if m is given by (jB.lSp . then (|B.10[) implies 

(B.18) TO± == (1 - co)mo + C±{fiw_) + C±{fiw+) 

= (1 - co)mo + C^ifJ.) ± fJ.w± 

= (1 - co)mo - (I - C„,)m + M^±- 
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From this we conclude that /i — m,±b^^ solves (|B.16|) . 
Conversely, if p, solves (|B.17|) , then set 

m{z)^mo{z) + — / /i(s)(w+(s) + w_(s))ri^(s,z), 

and the same calculation as in (JB.ISP implies rh± — fj.b±, which shows that m = 
(1 + co)^^m solves the Riemann-Hilbert problem (|B.1[) . D 

Note that in the special case 7 = we have mo{z) = (l l) and we can choose 
C as we please, say C = 00 such that cq = cq = in the above theorem. 

Hence we have a formula for the solution of our Riemann-Hilbert problem m{z) 
in terms of (I— Cu,)~^too and this clearly raises the question of bounded invertibility 
of I — Cw This follows from Fredholm theory (cf. e.g. [44]): 

Lemma B.4. Assume Hypothesis \B.l[ The operator I — Cw is Fredholm of index 
zero, 

(B.19) ind(I-C^)=0. 

Proof. Since one can easily check 

(B.20) (I - C^){1 - C_^) = (I - C_„)(I - C^) - I - T^, 

where 

Tw = T++ +T+_ +r-+ +T__, T^ta^U) = C^JCtral/W-^JW-^J, 

it suffices to check that the operators Taia2 are compact (34, Thm. 1.4.3]). By 
Mergelyan's theorem we can approximate 'w± by rational functions and, since the 
norm limit of compact operators is compact, we can assume without loss that w± 
have an analytic extension to a neighborhood of E. 

Indeed, suppose /„ G i^(S) converges weakly to zero. Without loss we can 
assume /„ to be continuous. We will show that ||Tu,/„||^2 — > 0. 

Using the analyticity of w in a neighborhood of S and the definition of C± , we 
can slightly deform the contour E to some contour E± close to E, on the left, and 
have, by Cauchy's theorem, 

T+ + fn{z) =^y {C{fnW-){s)w^{s))nds,z)- 

Now {C {fnW -)w -){z) ^ as n ^ 00. Also 

|(C(/„W-)w-)(z)| < co7isi||/„||i2||w_||Loo < const 

and thus, by the dominated convergence theorem, ||T_|__|_/„||2,2 — > as desired. 

Moreover, considering I — eCu, = I — Cew for < e < 1 we obtain ind(I — C^) = 
ind(I) = from homotopy invariance of the index. D 

By the Fredholm alternative, it follows that to show the bounded invertibility of 
1— Cw we only need to show that ker(I — Cw) = 0. The latter being equivalent to 
unique solvability of the corresponding vanishing Riemann-Hilbert problem in the 
case 7 = (where we can choose C, = 00 such that cq = cq = 0). 

Corollary B.5. Assume Huvothesis \B.l[ A unique solution of the Riemann-Hilbert 
problem (jB.ip with 7 = exists if and only if the corresponding vanishing Riemann- 
Hilbert problem, where the normalization condition is replaced by m{0) = (O '712), 
with TO2 arbitrary, has at most one solution. 
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We are interested in comparing a Riemann-Hilbert problem for which |jw||oo is 
small with the one-soliton problem, where 

(B.21) ||u;||oo = ||w+||l~(s) + ||w-||loo(s). 

For such a situation we have the following result: 

Theorem B.6. Fix a contour E and choose C; 7 = 7*; ^^ depending on some 
parameter t G R such that Hypothesis \B.1\ holds. 
Assume that w* satisfies 

(B.22) \\w*\\ao < Pit) 

for some function p{t) — > as t —> oo. Then (I — C'^t)^^ : i^(S) — > -^s(^) exists 



for sufficiently large t and the solution m{z) of the Riemann-Hilbert problems (jB.ip 
differs from the one-soliton solution TOq(z) only by 0{p{t)), where the error term 
depends on the distance of z to S U {C^^}- 

Proof. By boundedness of the Cauchy transform, one has 

||C„t|| < const 1 1 It;* I loo- 

Thus, by the Neumann series, we infcr that (I— Cj^jt)^^ exists for sufBcicntly large 
t and 

\\{I-C.^.)-^-I\\=0{p{t)). 
This imphes ||/i* — "JqIU^ — 0{p{t)) and Cq = 0{p{t)) (note jll = fi^ — mp). 
Consequently Cq — 0{p{t)) and ||/i* — m-olU^ — 0{p{t)) and thus m*-{z) — Wo(^) — 
0{p{t)) uniformly in z as long as it stays a positive distance away from S U {C^^}- 

D 
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